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Abstract 

We report the observation of B^ decays to the K^tt^tt^ final state using a data 
sample of 78 fb^^ collected by the Belle detector at the KEKB e+e^ collider. With 
no assumptions about intermediate states in the decay, the branching fraction is 
measured to be {SG.G'^^'^ it 3.0) x 10~^. We also search for B decays to intermedi- 
ate two-body states with the same K^tt^tt^ final state. Significant B signals are 



observed in the p{770) K^ and i^*(892)"'"7r channels, with branching fractions of 
(15.ll|3ti;5l|?) X 10-6 and (14.81^4-1:0-0:9) x 10"^ respectively. The first error 
is statistical, the second is systematic and the third is due to the largest possible 
interference. Contributions from other possible two-body states will be discussed. 
No CP asymmetry is found in the inclusive K^tt~tt^ or p^K'^ modes, and we 
set 90% confidence level bounds on the asymmetry of —0.12 < Acp < 0.26 and 
-0.18 < Acp < 0.64, respectively. 

PACS: 13.25.Hw, 14.40.Nd 



1 Introduction 



Recently, observations of large branching fractions for three-body charmless 
hadronic decays of B mesons have been reported by the B factory experiments 
[1,2,3,4,5,6]. In the mesonic decays B -^ Kn^n^ and B -^ KK^K~ , a large 
fraction of the decays proceed through intermediate two-body decay processes, 
such as B+ -^ ir*(892)°7r+, K*{892f -^ K+tt- and B+ -^ (l)K+, 4> -^ K+R-. 
However, higher mass K^n^ , tt+tt^ and K^K^ states may contribute but are 
not clearly identified due to limited statistics. Moreover, the broad K^K^ 
mass spectrum above 1.5 GeV/c^ in B^ -^ K^K^K^ suggests a large non- 
resonant B^ -^ K^K^K~ contribution. In the baryonic decay B^ -^ ppK~^, 
a simple phase-space model fails to describe the pp mass spectrum, which 
may be explained by a baryonic form factor model [7] or by an additional, 
unknown resonance around 2 GeV/c^. These studies of three-body decays have 
provided new information on the mechanism of B meson decay, and provide 
opportunities to search for unknown B meson decays and to understand the 
interference between them. In this paper we report on a study of B meson 
decays to K'^tt'tt'^, independently of possible intermediate states. In addition, 
we also present results of a search for quasi-two-body intermediate states. 
Inclusion of charge conjugate modes is always implied in this letter unless 
otherwise specified. The results are obtained from data collected by the Belle 
detector [8] at the KEKB asymmetric e~^e~ storage ring [9]. The data sample 
corresponds to an integrated luminosity of 78 fb"^ and contains 85.0 million 
BB pairs at the T(4S') resonance. 
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2 Apparatus and Event Selection 



The Belle detector is a large-solid-angle general purpose spectrometer based on 
a 1.5 T superconducting solenoidal magnet. Charged tracks are reconstructed 
with a three layer double-sided silicon vertex detector (SVD) and a central 
drift chamber (CDC) that consists of 50 layers segmented into 6 axial and 5 
stereo superlayers. The CDC covers the polar angle range between 17° and 
150° in the laboratory frame and, together with the SVD, gives a transverse 
momentum resolution of {apJptY = (0.0019 p^)^ + (0.0030)^, where pt and a^^ 
are in GeV/c. Charged hadron identification is performed using a combination 
of three devices: an array of 1188 aerogel Cerenkov counters (ACC) covering 
the momentum range 1-4 GeV/c, a time-of-flight scintillation counter system 
(TOF) for track momenta below 1.5 GeV/c, and dE/dx information from the 
CDC for particles with low or high momenta. Situated between these devices 
and the solenoid coil is an electromagnetic calorimeter (ECL) consisting of 
8736 CsI(T^) crystals with a typical front-surface cross-section of 5.5 x 5.5 
cni^ and a depth of 16.2 Xq. The ECL provides a photon energy resolution of 
{oEJEf = 0.0132 + (0.0007/E)2 + {0.008/ E^/y, where E and a^ are in GeV. 
An instrumented iron flux return outside the solenoid coil is used for muon 
and Kl detection. A detailed description of the Belle detector can be found 
in Ref. [8]. 

Charged tracks are required to come from the collision point and have trans- 
verse momenta, pt, above 100 MeV/c. The accepted tracks are then refitted 
with their vertex position constrained to the run-averaged profile of B meson 
decay vertices in the transverse plane. Charged K and tt mesons are identified 
by combining information from the CDC {dE/dx), the TOF and the ACC 
to form a K{tt) likelihood Lk{L.„). Discrimination between kaons and pions 
is achieved through the likelihood ratio L/^/(L^ + Lk)- The performance of 
the charged hadron identification is studied using a kinematically selected high 
momentum D*^ data sample, where D*^ -^ D^ir^, D^ -^ K~ti^ . We measure 
the pion and kaon identification efficiencies and their fake rates as functions 
of track momentum. The typical kaon and pion identification efficiencies for 
1 GeV/c tracks are (87.9 ± 0.6)% and (89.4 ± 0.6)%, respectively. The rate 
for true pions to be misidentified as kaons is (9.0 ± 0.5)%, while the rate for 
true kaons to be misidentified as pions is (10.0 ± 0.6)%. Charged tracks which 
are positively identified as electrons and muons are rejected. Candidate neu- 
tral pions are selected by requiring the two-photon invariant mass to be in 
the mass window 0.118 GeV/c^ < M(77) < 0.150 GeV/c^, corresponding to 
±2.5o" mass resolution with momentum above 2 GeV/c. The momentum of 
each photon is then readjusted, constraining the mass of the photon pair to be 
the nominal 7r° mass. To reduce the background from soft photons, each pho- 
ton is required to have energy above 50 MeV and the minimum ti^ momentum 
is 200 MeV/c. 



Candidate B mesons are identified using the beam constrained mass, Mbc = 
\/-^beam " -^B' ^nd the energy difference, AE = Eb - -Ebcam, where -Ebeam is 
run-dependent and determined from B -^ D^*^tx events, and Pb and Eb are 
the momentum and energy of the B candidate in the T(4S') rest frame. The pa- 
rameterizations of the signal in Mbc and AE are determined by a GEAN-based 
Monte Carlo (MC) [10] simulation of non-resonant B^ -^ K^-k^-k^ decays, and 
various quasi-two-body decays to the K^-k^-k^ final state. The signal param- 
eterization is verified using the data and MC samples of 5^ -^ D^ir'^, D^ -^ 
K^TT~Ti^ candidates. Our MC overestimates the Mbc resolution by 8% but 
underestimates the AE resolution by 9% to 15%, depending on the kinemat- 
ics of the K^T\:~Tx^ events. The MC based signal probability density functions 
(PDF) are readjusted accordingly. 

The Gaussian width of the signal in Mbc is about 3.0 MeV/c^, which is pri- 
marily due to the beam energy spread. The AE distribution is found to be 
asymmetric, with a tail on the lower side due to 7 interactions with material 
in front of the calorimeter, and shower leakage out of the back side of the 
crystals. As a result, the AE resolution and the tail distribution strongly de- 
pend on the tt" energy; the AE width ranges from 20 MeV to 33 MeV. In the 
inclusive K^tt^tt^ study, since the tt^ energy distribution for the signal is not 
known a priori, the data is divided into three samples: -P(7r°) < 0.5 GeV/c, 0.5 
GeV/c < P(7r°) < 1.5 GeV/c and P(7r°) > 1.5 GeV/c. The AE distribution 
in each sample is modeled with a Crystal Ball lineshape [11] with parameters 
determined from MC. Events with Mbc > 5.2 GeV/c^ and \AE\ < 0.3 GeV 
are selected for the final analysis. The signal region is defined as Mbc > 5.27 
GeV/c^ and -0.10 GeV < AE < 0.08 GeV. Events located in the region 
Mbc < 5.265 GeV/c^ are defined as sideband events and are used for back- 
ground studies. When more than one B^ candidate is found in an event, the 
candidate having the smallest sum of the x^ from the vertex fit and 7r° mass 
constrained fit is selected. 

The dominant background for three-body B decay events comes from the 
e~^e~ — > qq continuum, where q = u,d,s or c. In order to reduce this back- 
ground, several shape variables are chosen to distinguish spherical BB events 
from jet-like continuum events. Five modified Fox- Wolfram moments [12] and 
a measure of the momentum transverse to the event thrust axis {S±) [13] 
are combined into a Fisher discriminant. The PDFs for this discriminant and 
cos6b, where 6b is the angle between the B fiight direction and the beam 
direction in the T(45') rest frame, are obtained using events in the signal and 
sideband regions from MC simulations for signal and qq background. These 
two variables are then combined to form a likelihood ratio TZ = Cs/{Cs + Cqq), 
where Cs(gg) is the product of signal {qq) probability densities. Continuum 
background is suppressed by requiring TZ > 0.9, based on a study of the signal 
significance {Ns/y^Ns + Nb) using a MC sample, where Ns and Nb are signal 
and background yields, respectively. This TZ requirement retains 45% of the 



signal and removes 97% of the continuum events. The effect of the TZ cut is 
studied by comparing B'^ — * D^7r~^ in data and MC, for different values of TZ. 
A systematic error of 3% is obtained for the TZ cut. 
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Fig. 1. Dalitz plot of B^ -^ K^tt^tt^ candidates from (a) the B signal region and (b) 
the Mbc sideband region. The enhancement around M'^{K~^7r~) = 3.5 (GeV/c^ 
from B^ -^ D^TT^ events. 
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Figure 1 shows the Dalitz plot distribution of K^tt^ti^ candidates in the AE — 
Mbc signal region and Mbc sideband region. K^-n-^n^ candidates populate 
the three edges of the Dalitz plot, indicating the existence of quasi-two-body 
intermediate states. Moreover, there is an enhancement near M'^{K 
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TT 
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(GeV/c ) in Fig. 1(a), which is due to the decay B" -^ D 
To restrict the study to charmless B decays, events with a K'^n 
50 MeV/c^ of the nominal D^ mass are rejected. 
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3 Inclusive K'^n 7r° Yield 



The final signal yields are obtained from fits to the AE and Mbc distributions. 
In addition to continuum background, the final sample contains background 
from T(4S') -^ BB events. The AE and Mbc shapes of this BB background 
are modeled with smooth histograms, generated from a large GEANT based 
BB MC sample, which includes b ^ c transitions and charmless B decays. 
The continuum AE background shape is modeled by either a first or second 
order polynomial, determined from the Mbc sideband data. The continuum 
Mbc background shape is modeled with an ARGUS function [14] with param- 
eters determined from events outside the AE signal region. One-dimensional 
binned likelihood fits to AE and M^c are performed using signal, continuum 
background and BB background PDFs for events in the Mbc and AE signal 
region, respectively. Since the Mbc shapes from BB background are difficult 



Table 1 

Fit results for inclusive K^tt^tt^ events. Column 3 and 4 list the signal yield and 
the statistical significance. The BB fraction is fixed to the MC expectation in the 
Mbc fit, while it is allowed to float in the AE fit. The last row shows the sum of the 
three AE fits. 



vr" momentum (GeV/c) 


Signal Yield 


Significance 


Mbc fit 


No Requirement 


369 ± 35 


11.4 




p< 0.5 


3llii±4 


2.6 


A^fit 


0.5 <p< 1.5 


93 ± 22l^ 


4.6 




p > 1.5 


262 ± 36 ± 13 


7.8 




Sum 


386 ± Utlt 


9.4 



to distinguish from signal shapes, the signal yields are estimated using the 
AE fit and cross checked by the Mbc fit, where the BB background fraction 
is fixed to the MC expectation. 



Table 1 summarizes the fit result of the inclusive K 
statistical significance (E) defined as 



'71 7r° 



sample with the 

21n(£o/'Cmax), where Cq and £rnax 
denote the likelihood values at zero yield and the best fit numbers, respec- 
tively. The sum of the signal yield from the AE fits to the three subsamples, 
386 ± 44, is consistent with the yield from the Mbc fit, 369 ± 35, which has a 
smaller signal efficiency than the AE fit due to the tighter AE requirement. 
The corresponding projections of the fits are shown in Fig. 2. Furthermore, a 
consistent result is obtained when the BB fraction is fixed according to the 
MC expectation. 




0.0 
AE (GeV) 

Fig. 2. (a) Mbc and (b) AE distributions for inclusive K^ir^i:^ events. Solid lines 
in the figures represent the result of the fit and the signal contribution. Dashed lines 
show the total background contribution, while dotted lines indicate the continuum 
contribution. The AE lines in (b) show the sum of the fit results to the three 
subsamples. A sizeable 5° -^ K+vr" and _B+ -^ i^+7r° feed-down at AE > 0.2 GeV 
is found in the MC. This contribution is included in the fit. 
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Fig. 3. -B yields from AE fits as a function of (a) K^tt^, (b) K^tt^ and (c) vr^Tr". 
Each two-body mass is examined after requiring the other two-body masses to 
be large (> 1.6 GeV/c^ for Ki: and > 1.1 GeV/c^ for vr^vr'^). Dotted points are 
data; the superimposed curves in (b) and (c) are the projection curves based on 
i^*(892)"'" and /9(770)~, The enhancements between 1.1 to 1.5 GeV/c^ in (a) and 
(b) are modeled with Breit-Wigner functions. 

4 Two-body Intermediate States 



We perform a search for quasi-two-body decays in the K'^tt'tt^ final state, 
including B decays to a pseudoscalar {K or tt) and a vector meson (p(770) 
or K*{892)) and other possible intermediate states with higher mass reso- 
nances. Three pseudoscalar- vector (PV) modes are considered: i^*(892)°7r*', 
i^*(892)"'"7r~, and p{770)^ K^ . Figure 3 shows the B signal yields from the 



AE fit as functions of K^ 



7r 



K^TT^ and vr vr*' masses. To eliminate cross- 



talk between decay modes, each two-body mass is examined after requiring 
the other two-body masses to be large (> 1.6 GeV/c^ for Kn and > 1.1 
GeV/c^ for 7r~7r''). The AE signal PDFs are obtained from MC simulations 
of 5° -^ ir*(1430)°7r°, B^ -^ K*{U30)+7r- and B^ -^ p{770)-K+ for the 



K^Ti , K~^7T^ and tt 7r° cases, respectively. 

In the K^Ti^ sample, a large enhancement is observed between 1.0 and 1.6 
GcV/c^, peaking around 1.2 to 1.4 GeV/c^. More structure is observed in the 
K^TT^ spectrum. An enhancement is seen in the i^*(892)+ mass region, in the 
region from 1.2 to 1.4 GeV/c^ and possibly between 1.8 and 2.1 GeV/c^. In 
the 7r~7r° sample, a clear excess is seen in the p(770)~ signal region. Although 
the enhancement between 1.1 and 1.6 GeV/c^ is observed in both K^tt^ and 
K^-K^ spectra, these higher mass Kn states cannot be identified without per- 
forming an angular analysis that requires much more data. Possible candidates 
are ir*(1410), ir*(1430) and ir*(1430). Earlier studies of B+ -^ K+tt+tt- [1] 
and B^ — > K^7r~^7i~ [5] decays also observed large quasi-two-body B decays 
with K*{Kir) mesons in the final state. If the enhancements observed in the 
i^+TT^ mass spectrum are due to such K* mesons, one would expect the same 
resonances to appear in the K^7i~^ mode. The first two enhancements seen in 
the K'^ir'^ mode (Fig. 3(b)) are indeed observed in the B^ -^ K^ti^ti^ anal- 
ysis [2]. However, the third enhancement, between 1.8 and 2.1 GeV/c^, does 
not appear in the K'^7r~^7r~ mode. This enhancement, which has a signal yield 
of 22^*^8 events and a significance of 3a, may be either a statistical fiuctuation, 
or originate from i^2(1820), 7^4(2045) or from a doubly Cabibbo suppressed 
D^ decay. More data are needed to clarify the current situation. 

To further understand the possible resonances, we study the distributions 
of cos Oh, where the helicity angle 6h is defined as the angle between the 
direction of the candidate B meson and the K~^ (tt^) direction in the K*{p) 
rest frame. Figure 4 shows the B yields as a function of cos6h for events 
in the K*{892)~^, p{770)~,K*{Kn) signal region. The cos 9h distributions for 
the first two modes are consistent with those of B^ -^ pseudoscalar vector 
(PV) decays, as expected for 5° -^ fs:*(892)+7r' and 5° -^ p{770yK+. Note 
that the asymmetry in the cos6h distributions is due to the inefficiency of 
low momentum 7r° reconstruction. Since tt^s from B^ -^ p~ K^ decays are 
more energetic than those from B^ -^ K*(892)+7r~, the asymmetric effect is 
less pronounced. The cos Oh distributions for the K^ and K*"*" modes favor a 
scalar behavior. 



5 Yield for Various States 



We measure B^ -^ K^-k^-k^ decay rates for the three PV modes, events in 
the two -ft'*vr regions, and the central region of the Dalitz plot with two- 
body masses above 2.0 GeV/c^. Candidate i^*(892) and p(770)~ mesons are 
identified by requiring the K^n'^'^^ and 'k~ti'^ masses to be in the range 820- 
980 MeV/c^ and 570-970 MeV/c^, respectively. To further reduce background, 
a selection of | cos6h\ > 0.3 is applied to the vector meson candidates. In the 
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Fig. 4. i? signal yields as a function of cos 6h for events in (a) i^*(892)+, (b) p(770) , 
(c) K*^ and (d) i^^"*" regions. Points are data and histograms are the expectations 
from (a) B° -^ K*{892)+tt-, (b) B° -^ p{770)- K+ and (c,d) B decays to a scalar 
and a pseudoscalar meson. 

two Kf^+h regions, we require 1.1 GeV/c^ < M(ir+7r-(o)) < 1.6 GeV/c^. S 
meson candidates are then selected from the inclusive K^ti^ti^ events after 
applying all analysis cuts, including the appropriate two-body mass vetos to 
avoid cross talk. The signal PDFs are obtained from MC simulations for all 
six channels, where a scalar hypothesis is used to model K*tt. 



Figure 5 shows the Mbc and AE distributions, and their corresponding fit 
curves, for the three PV modes. No signal yield is seen in the K*^n^ chan- 
nel but significant signals are observed for the K*~^7i~ and p~K^ modes; the 
yields measured from the AE fits are 38 ±11 and TTl^y events, with statistical 
significances of 3.8(T and 4.9cr, respectively. As for events in the central region 
of the Dalitz plot and the two K^ regions. Fig. 6 shows their Mbc and AE 
distributions with the fit curves superimposed. Based on the AE fit, there are 
67 ± 17 and 52 ± 15 signal events in K^^ti^ and i^*"'"7r~, respectively. Since 
events in these K^n regions cannot be positively identified, their reconstruc- 
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Fig. 5. Mbc and AE distributions for (a,b) K*°Tr°, (c,d) K*+tt- and (e,f) p+A'~ 
events. The superimposed solid curves represent the corresponding fits. Dashed 
lines are the background projections and the dotted lines show the continuum qq 
contribution. Enhancements in the projection curves at AE > 0.2 GeV in (d) and 
(f) are due to B ^ Ktt decays. 

tion efficiencies are determined without assumptions about the intermediate 
two-body states. Although a yield of around 20 events is obtained from both 
the AE and Mbc fits for the central region of the Dalitz plot, the statistical 
significance is below 3a and, hence, an upper limit is reported. The recon- 
struction efficiency is obtained from a phase-space decay model. 

The number of feed-across events from high mass K* to the i^*(892) region is 
estimated from the B yields in the 1.1 GeV/c^ < M(A'+7r°) < 1.6 GeV/c^ re- 
gion, and the Breit-Wigner distribution, modeled with a mass of 1.326 GeV/c^ 
and a mean of 252 MeV/c^. We find a contribution of 1 event. Assuming no 
interference, the K*{892)~^7t~ yield is estimated to be 37 ± 11. The possible 
effect of interference is studied using a Monte Carlo simulation which assumes 
the three PV decays and two K*7i states. We compare the yields with and 
without interference. After varying the relative phase of each channel, and tak- 
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Fig. 6. Mbc and AE distributions for events in different regions of the Dalitz plots: 
(a,b) the central region, (c,d) K*^tt^ region and (e,f) K*'^7r~ region. The superim- 
posed solid curves represent the corresponding fits. Dashed lines are the background 
projections and the dotted lines show the continuum qq contribution. Enhancements 
in the projection curves at AE > 0.2 GeV in (d) and (f) are due to i? — > Kir decays 



ing into account the non-uniform reconstruction efficiency over the Dalitz plot, 
the largest deviation is ^1%° for K*{892)+tt- and ~!:lf% for p-K+. These two 
numbers are used to estimate the systematic error arising from interference. 



6 Systematic Uncertainties 



The systematic error for each signal yield is estimated by varying each param- 
eter of the ffi functions by ±lo" from the measured values. The shifts in signal 
yield are then added in quadrature. The typical fit systematic error for the 
inclusive decay is around 4%. Signal efficiencies are ffist obtained from MC 
simulations, and then corrected by comparing data and MC predictions for 
other processes. The efficiency for the inclusive K^tt^ti^ signals is estimated 
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Table 2 

Summary of the B^ -^ K^tt^tt^ search. We present signal yields, efficiencies, and 
their statistical significances for the inclusive mode, the three intermediate chan- 
nels of i? ^ PV decays and the other three regions of the Dalitz plot. Branching 
fractions and/or upper limits are shown in the last two columns. In the branching 
fractions, the first error is statistical and the second systematic. The third error for 
i^*(892)+7r~ and p~K^ corresponds to the largest uncertainty from the interference 
between different states. The last channel, i^+7r~7r^^, indicates the non-resonant 
S° -^ K+'K~'K^ decay. 



Channel 


Yield 


Eff.(%) 


Sig. 


BF (10-6) 


UL (10^6) 


K+vr-vrO 


386 ± 44 


12.4 


9.4 


36.6l^-^ ± 3.0 


- 


K*{smfTi^ 


9 + 11 
^-10 


6.7 


0.3 


0.41};? ±0.1 


3.5 


K*(892)+7r- 


37 ±11 


2.9 


3.8 


1/1 0+4.6+1.5+2.4 
l"t-0_4.4-l. 0-0.9 


- 


p-K+ 


77+18 
''-17 


6.0 


4.9 


ic; 1+3.4+1.4+2.0 
I'^-^-S. 3-1. 5-2.1 


- 


KflT^ 


67 ±17 


12.9 


4.2 


f. 1+1.6+0.5 
•3-1-1.5-0.6 


- 


K*+TT- 


52 ±15 


11.9 


3.7 


5.1 ± 1.5l|]-6 


- 


K^t^'Ar 


22+10 


4.1 


2.5 


r 7+2.7+0.5 
^•'-2.5-0.4 


9.4 



from the weighted sum of the efficiencies for the possible two-body interme- 
diate states shown in Fig. 3, where the sub-decay branching fraction for each 
two-body state is not included. The uncertainty on this inclusive efficiency 
is 5%, determined by checking the reconstruction efficiencies on various two- 
body modes. The tt^ reconstruction efficiency is verified by comparing the vr^ 
decay angular distribution with the MC prediction, and by measuring the ratio 
of the branching fractions of two rj decay channels: r] ^ 77 and rj -^ tt^ti'^ti^. 
The typical systematic error for n^ detection is 3%. The systematic errors on 
the charged track reconstruction are estimated to be ~ 2% using partially 
reconstructed D* events, and verified by comparing the ratio of r^ ^ Tr+Tr^TT^ 
to ?7 ^ 77 in data with MC expectations. The final systematic errors on the 
reconstruction efficiencies, including charged particle and vr'^ detection, parti- 
cle identification and the 7^ cut, range from 7 to 14% for the quasi-two-body 
decays and is 8.2% for the inclusive i^+7r~7r° channel. 



7 Branching Fractions 



Table 2 summarizes the ffi results for each reconstructed decay channel. The 
branching fractions and upper limits are calculated assuming that B^B~ 
and B^B^ are produced with equal probability. The systematic errors on the 
branching fractions combine the systematic errors for the reconstruction effi- 
ciencies and the AE fit with the uncertainty from the number of BB events. 
Since no signal is seen in the i^*°7r° mode, and the signal yield in the central 
region of the Dalitz plot is not significant, upper limits are computed at the 
90% confidence level (C.L.) based on the observed number of events in the sig- 
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Fig. 7. AE distributions for the (a) K+7r"7r°, (b) K"7r+7r°, (c) p-K+ and (d) 
p^K~ samples. The superimposed sohd curves represent the corresponding fits and 
signal projections. Dashed lines are the background projections and dotted lines 
show the continuum qq contribution. Enhancements in the projection curves at 
AE > 0.2 GeV are due to i? ^ Kn decays 



nal region, and the background level found by the fit; both the statistical and 
systematic errors are taken into account [15]. With 85.0 million BB events, we 
measure the branching fraction to be (36.6t4;i(stat.)±3.0(syst.)) x 10~^ for the 
inclusive B^ — > K~^TT~7r^ decay, without assumptions about intermediate two- 



body states. The branching fractions of B -^ K*(892)' 



TX 



and 5° -^ p-K+ 

decay are measured to be {U.Sttitlitli) x lO^^ and il5.ltli1:\itli) x 10"^ 
where the first error is statistical, the second is systematic and the third cor- 
responds to the largest uncertainty from the interference between different 
states. Finally, the B^ -^ K~^tt~tt'^ decay branching fractions in the K*'^7i^ and 
K*+TT- regions are measured to be (6.ll||t[!:^) x lO^^ and (5.1±1.5+g:^) x 10"^ 
respectively. 
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8 Search for CP Violation 



Using the large signals observed in the inclusive B^ -^ K^ti^ti^ and B^ -^ 
p~ K^ modes, we search for direct CP violation by dividing the data into 
two subsets, according to the charge of the kaon. The asymmetry, defined as 
^cp = jv^^Jv^' ^^ then computed using the B signal yields obtained from 
A£^ fits. Following the same fitting procedure, we observe 179^30 K^n^Ti^ 
events and 2071*^3^ K~n^Ti^ events, while the /S.E fit yields 30lj;^ p~ K^ events 
and 47]!^|2 P^K~ events (see Fig. 7). The possible reconstruction bias in Acp is 
studied by checking the inclusive D^ — * K'^tt" and D^ -^ K^-k^ yields in data. 
The obtained systematic error is 0.5%. Adding this 0.5% error in quadrature 
with the fitting systematic error, obtained by varying each parameter in the 
PDFs by lo", gives the total systematic error. Finally, the CP asymmetry is 
calculated to be Acp = 0.07 ± 0.11 ± 0.01 for the inclusive mode, and Acp = 
0.22+[]:|+n^ for B^{W) -^ p^K^. We also set 90% confidence intervals on 
the asymmetry of —0.12 < Acp < 0.26 for the inclusive mode, and —0.18 < 
Acp < 0.64 for 5° -^ p^K^. 



9 Conclusions 



In summary, we have studied the charmless hadronic decays, B^ -^ K'^tt 7r°, 
which is observed for the first time. Our results show that the branching 
fraction of B^ -^ K+n-n^ is (64 ± 10)% and (78 ± 15)% of that of B+ -^ 
K+TT-TT+ [2,4] and fi° -^ K°tt-tt+ [2,5], respectively. The K+n-n^ signal 
candidates populate the edge of the Dalitz plot, indicating the existence of 
quasi- two-body states. For the K^ti~'k'^ final state, we have observed signals 
in the K*^tt~ and p~K^ samples but no significant K*^ti^ signal is seen. The 
p^K^ branching fraction is close to the K*^7t^ branching fraction, where 
our measurement is consistent with the earlier CLEO result [5]. However, our 
p~ K^ result is twice that of BaBar's measurement [6]. We also report the 
B decay rates in other regions of the Dalitz plot without assumptions about 
the presence of two-body intermediate states. In the future, significantly more 
data will be collected at Belle, which will enable us to perform a full Dalitz 
analysis, allowing us to identify other quasi-two-body states and extract their 
relative phases. Finally, we performed a search for direct CP violation in the 
inclusive and B^ -^ p~K^ channels. No evidence of CP violating asymmetry 
is seen and 90% C.L. limits on Acp are set. 

Acknowledgements 

We wish to thank the KEKB accelerator group for the excellent operation 
of the KEKB accelerator. We acknowledge support from the Ministry of Ed- 



16 



ucation, Culture, Sports, Science, and Technology of Japan and the Japan 
Society for the Promotion of Science; the Australian Research Council and 
the Australian Department of Education, Science and Training; the National 
Science Foundation of China under contract No. 10175071; the Department 
of Science and Technology of India; the BK21 program of the Ministry of 
Education of Korea and the CHEP SRC program of the Korea Science and 
Engineering Foundation; the Polish State Committee for Scientific Research 
under contract No. 2P03B 01324; the Ministry of Science and Technology 
of the Russian Federation; the Ministry of Education, Science and Sport of 
the Republic of Slovenia; the National Science Council and the Ministry of 
Education of Taiwan; and the U.S. Department of Energy. 



References 

[I] Belle Collaboration, A. Garmash et al, Phys. Rev. D 65 (2002) 092005. 

[2] Belle Collaboration, A. Garmash et a/., Phys. Rev. D 69 (2004) 012001. 

[3] Belle Collaboration, M. Wang et al, Phys. Rev. Lett. 92 (2004) 131801. 

[4] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 91 (2003) 051801. 

[5] CLEO Collaboration, E. Eckhart et al, Phys. Rev. Lett. 89 (2002) 251801 . 

[6] BaBar Collaboration, B. Aubert et al, Phys. Rev. Lett. 91 (2003) 201802. 

[7] C.-K. Chua, W.-S. Hon, and S.-Y. Tsai, Phys. Rev. D 65 (2002) 034003; H.-Y. 
Cheng and K.-C. Yang, Phys. Rev. D 66 (2002) 014020. 

[8] Belle Collaboration, A. Abashian et al, Nucl. Inst, and Meth. A 479 (2002) 
117. 

[9] S. Kurokawa and E. Kikutani, Nucl. Inst. Meth. A 499 ( 2003) 1, and other 
papers included in this Volume. 

[10] R. Brun et al, GEANT 3.21, CERN Report No. DD/EE/84-1 (1987). 

[II] Crystal Ball Collaboration, J.E. Gaiser et al, Phys. Rev. D 34 (1986) 711. 

[12] The Fox- Wolfram moments were introduced in G. Fox and S. Wolfram, Phys. 
Rev. Lett. 41 (1978) 1581. The modified Fox- Wolfram moments were described 
in K. Abe et al, Phys. Lett. B 511 (2001) 151. 

[13] CLEO Collaboration, R. Ammar et al, Phys. Rev. Lett. 71 (1993) 674. 

[14] ARGUS Collaboration, H. Albrecht et al, Phys. Lett. B 241 (1990) 278. 

[15] For the i^*(892)''7r'' limit we follow the approach of G.J. Feldman and 
R.D. Cousins, Phys. Rev. D 57 (1998) 3873, using the implementation of 
J. Conrad et al, Phys. Rev. D 67 (2003) 012002 to include the systematic 



17 



errors. For the i^'*~7r~7rjy^ mode this method returns a two-sided interval 
[2.0, 9.4] X 10"^. Because the significance is less than three, we discard the lower 
edge of this interval and (conservatively) quote an upper limit of 9.4 x 10~^. 



